Quantum phase transition in a two-channel-Kondo quantum dot device 
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We develop a theory of electron transport in a double quantum dot device recently proposed 
in Ref. Q for the observation of the two-channel Kondo effect. Our theory provides a strategy for 
tuning the device to the non-Fermi-liquid fixed point, which is a quantum critical point in the space 
of device parameters. We explore the corresponding quantum phase transition, and make explicit 
predictions for behavior of the differential conductance in the vicinity of the quantum critical point. 
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I. INTRODUCTION 

The magnetic screening of a localized spin by spins 
of itinerant electrons^ leads to the Kondo effect - an 
anomaly in low-temperature conduction properties. This 
screening becomes effective below some characteristic 
temperature, the Kondo temperature Tk- Above Tk 
electrons are weakly scattered by the magnetic impurity, 
but below Tk the scattering becomes strong. In the sim- 
plest Kondo systems, only one electron mode (the s-wave 
mode, say) participates in the screening of a localized 
spin with S = 1/2. In this case, the low-temperature 
electronic properties are adequately described by Fermi 
liquid theory, and the thermodynamic and transport 
characteristics are analytical functions of T /Tk- In more 
complicated systems (such as, e.g., paramagnetic metals) 
many electron modes may participate in screening of an 
S = 1/2 localized moment 4 . The peculiarities of such 
a "multichannel" Kondo model were long recognized 4 *^. 
At the same time it was understood that even a small 
deviation from symmetry between channels leads at low 
temperatures to the Kondo screening by just one channel, 
the one for which the exchange integral with the impurity 
is the largest. 

The peculiarity of a symmetric multichannel Kondo 
problem is in its non-Fermi-liquid (NFL) behavior at 
low temperatures^. The low-temperature asymptotes of 
the thermodynamic and transport characteristics display 
power-law behavior with fractional values of the expo- 
nents. A complete temperature dependence of the ther- 
modynamic characteristics (such as the local spin sus- 
ceptibility) is known now from the exact Bethe-ansatz 
solution of the Kondo problem*^. Details of the low- 
temperature electron scattering problem were also un- 
derstood in the framework of conformal field theor^S. 

Experimental observation of the non-Fermi-liquid be- 
havior in a Kondo system, however, is difficult because 
the channel symmetry is not "protected" - in general, 
there are no conservation laws prescribing such a symme- 
try. This has lead to various propositions to observe such 
a behavior in systems where the role of spin is taken over 
by another degree of freedom, while the "real" spin la- 



bels the channels, making the channel symmetry robust. 
One such idea deals with an atomic defect which occupies 
two equivalent lattice sites, thus forming a pseudospin 10 . 
However, the equivalence of sites is not a protected sym- 
metry; its violationii, equivalent to a "Zeeman splitting" 
of the pseudospin states, destroys the Kondo effect. 

Another object which under certain conditions can be 
described by the two-channel Kondo model (2CK) model, 
is a large quantum dot, or a metallic island connected by 
a single-mode channel to a conducting electrode^. If 
one neglects the finite level spacing in the island, then 
a pseudospin labeling of the charge states of the island 
may be introduced, while real spin again plays the part of 
the channel index. In this setup the degeneracy with re- 
spect to the pseudospin orientation is easily achieved by 
tuning the gate voltage to the vicinity of the Coulomb 
blockade degeneracy point. At temperatures T higher 
than the level spacing SE in the island, the system is 
then described by the 2CK modeli^. Since Tk for this 
system can be of the order— of the charging energy Eq , 
while typically SE -C Ec, the NFL regime is easily re- 
alized. When an additional electrode is attached to the 
island, one can study the transport properties of the re- 
sulting device. The disadvantage of such realization of a 
2CK system is that there is no mapping between the con- 
ductance across the island*^ and the electron scattering 
cross-section in the generic two-channel Kondo model^. 

Small quantum dots with large level spacing have 
proved to be suitable for the observation of the Kondo 
effect. In the usual geometry consisting of a dot with 
two attached electrodes, however, only the conventional 
Fermi-liquid (FL) behavior is observable at low temper- 
atures. The reason lies in the structure of the matrix of 
exchange constants that couple the dot's spin to the spins 
of itinerant electronsi§*i£. Typically, the eigenvalues of 
this matrix are vastly differen1pi&, and their ratio is not 
tunable by conventional means. 

A device that circumvents this problem was proposed 
recently in Ref. 01 and involves several dots. A two-dot 
device is sufficient for the realization of the 2CK model. 
The key idea of Ref. [j] is to replace one of the electrodes 
in the standard configuration by a very large quantum 
dot 2, see Fig. ^ characterized by a level spacing 8E2 
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and a charging energy E^. At T 3> SE2, particle-hole 
excitations within this dot are allowed, and electrons of 
dot 2 participate in the screening of the smaller dot's 
spin. At the same time, as long as T <C E2, the number 
of electrons in the dot 2 is fixed. As a result, the electrons 
in dot 2 provide for a separate channel which does not 
mix with the channels provided by the electrodes L and 
R. In this case, the exchange constants for two channels 
may be tuned to become equai: the asymmetry between 
the channels is controlled by the ratio of the conductances 
of the dot-leads and dot-dot junctions. 




FIG. 1: Device proposed in Ref.Q. Level spacing in the larger 
dot (2) must be negligibly small to allow for the NFL behavior 
of the device at low temperatures. 

In principle, a setup having just one lead and two 
dots would allow one to study thermodynamic properties, 
such as magnetic susceptibility, in the 2CK regime. The 
existing technology 1 -, however, enables one to measure 
transport rather than thermodynamic properties. There- 
fore, two leads are needed to perform conductance mea- 
surements. In this paper, we assume that one of the elec- 
trodes is coupled weakly to the small dot and serves as a 
probe of the 2CK system formed by the two dots and the 
remaining electrode. We propose a detailed strategy for 
tuning the device to the NFL regime, and discuss various 
manifestations of NFL-related physics in the transport 
properties of the system. 



II. THE MODEL 

According to the discussion above, the device we con- 
sider consists of two quantum dots coupled to two con- 
ducting leads via single-mode junctions. The model 
Hamiltonian of such a device can be written as a sum 
of three parts, 



H = H d + H l +H t 



(1) 



The first term here, Hd, describes an isolated system of 
two quantum dots, 1 and 2, connected via a single mode 
junction, 



y 2ks 



ks \ ks / 

EMfcA+ILc.) ( 2 ) 



The last two terms in Eq. £[} represent the free electrons 
with spin s = ±1 in leads R and L, and the tunneling 
between the leads and dot 1, see Fig.^ 



H * = XX c L s c afc s > a = R,L; (3) 

aks 

Ht = E^LA+H-c. (4) 

aks 

In Eq. the smaller dot (dot 1) is described by a single- 
level system equivalent to the Anderson impurity model. 
The parameter E\ represents charging energy, while the 
parameter N is adjustable by tuning the potential on 
the capacitively coupled gate electrode. We neglect the 
finite level spacing SE2 in the dot 2, but account for its 
finite charging energy E2 (we do not write explicitly the 
gate potential applied to the dot 2, as it corresponds to 
a trivial shift of the chemical potential). 

Since the relevant energies (ui < Tk) for the Kondo 
effect are negligibly small compared to the Fermi energy, 
the electronic dispersion relation £/. in Eqs. ©,£3 can be 
linearized: = vpk, where k is measured from the Fermi 
momentum fejr. The linearization leads to an energy- 
independent density of states v, which will be assumed 
throughout this paper. Finally, we treat the tunneling 
amplitudes t2,tn,tL as real numbers and neglect their 
dependences on k. This is well justified for relevant values 
of k, \k\ < T/v F . 

Instead of working with the operators CR t L, it is con- 
venient to introduce their linear combinations ipo,i, 



^Oks 



cos sin #o 
— sin #0 cos 8q 



CRks 
Chks 



(5) 



where the angle #0 is determined by the equation 

tan 6»o = t L /t R . (6) 

(So far there are no restrictions on the value of t^/tR.) 
The Hamiltonian then assumes the "block- 

diagonal" form 



H 

Ho 

Hi 



H {il>o} + Hi {tp u tp 2 ,d}, 

E ^k^Oks^Oks' 
ks 

H d {^2,d} + J2^k4ksAks 

ks 

ks 



(7) 
(8) 

(9) 



where Hd{il>2,d} is given by Eq. J5J, and t\ = \Ji 2 L + i l R . 

At low energies (T -C £4,2) the Hamiltonian Hi in- 
volving the i/ji and 1^2 operators, see Eq. |J3J, can be 
simplified further. Indeed, at N ps 1 the small dot is 
occupied by a single electron, and, therefore, carries a 
spin 5 = 1/2. The tunneling terms in Eqs. @ and © 
mix the states with a single electron in in dot 1 with 
states having or 2 electrons in that dot. Because of the 
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high energy cost (~ E\), these transitions are virtual, 
and, provided that the conductances of the correspond- 
ing junctions are small, can be taken into account pertur- 
batively in the second order in tunneling amplitudes. A 
newi and important element here compared to the con- 
ventional treatment of the Anderson impurity model is 
that at T <C E 2 only those excitations that conserve the 
number of electrons in dot 2 are allowed. The result- 
ing effective Hamiltonian which acts within the strip of 
energies \u>\ < mm{Ei, E 2 }, has the form of the 2CK 
model 5 . 6 . 7 . 8 . 9 . 10 , 

H 2C k = £ ^\ k J lks + £j 7 ( s 7 ■ S) + BS Z . (10) 

7/cs 7 

Here the channel index 7 = 1 and 7 = 2 represents the 
leads and dot 2, respectively, S is the spin-l/2-operator 
describing the doubly-degenerate ground state of dot 1, 

kk' ss' 

is the spin density in channel 7, and cr = (a x , a v , a z ) 
are the Pauli matrices. The exchange amplitudes J 7 in 
Eq. (|1(J|> are estimated as 

vJ 1 =Avt*/E 1 . (11) 

In derivation of Eq. (|l(Jfl . we assumed that the gate 
voltage is tuned precisely to N = 1 (which corresponds 
to a particle- hole symmetric situation). As we discuss 
in Section [V] below, this assumption does not lead to 
qualitative changes in the results. We also included in 
the Hamiltonian the effect of an external magnetic field 
(hereinafter we omit the Bohr magneton /i^; the field B 
is measured in the units of energy). 

III. TUNNELING CONDUCTANCE 

In order to study the out-of-equilibrium transport 
across the device we add to our Hamiltonian a term 

H V = ^ (Nl - Nr) , #« = £4*.c ak „ (12) 

ks 

which describes a finite bias voltage V applied between 
the left (a = L) and right (a = R) electrodes. The dif- 
ferential conductance dl/dV can be evaluated in a closed 
form for arbitrary V when one of the leads, say L, serves 
as a weakly coupled probei 7 -, i.e., tz -C £r. Under this 
condition the angle do in Eqs. J5J and © is small: 

0o » t L /t R « 1. (13) 

Application of the transformation Eq. JSJ to Eq. (|12H . 

yields, to the linear order in 9q, 

Hv = ^ (No- N 1 )+eVe J2(4 ks ^i ks +^ ( 14 ) 

ks 



where 

ks ks 

The first term on the right-hand side of Eq. i|14[l can be 
interpreted as a voltage bias between the reservoirs of Gl- 
and 1-particles, cf. Eq. i|12f> . while the second term has 
an appearance of the fc-conserving tunneling. Since the 
tunneling amplitude is proportional to the small parame- 
ter #0^1; scc Eq. ((13)1 , one can use perturbation theory 
to calculate the current across the devicei 7 -. 

Similar to the representation of Hy in the form of 
Eq. i|14|) . the current operator 

/= if(**-**) 

also splits naturally into two contributions, 

I = J + SI. (15) 

Here 

/o = H(#i-*o) (16) 

ks 

is a current between the reservoirs of 0- and 1-particles, 
and 

It is easy to show^i that in the leading (second) order in 
9q the operator SI does not contribute to the average cur- 
rent across the device. The remaining contribution (Jo) 
corresponds to the ^-conserving tunneling between two 
bulk reservoirs containing 0- and 1-particles, see Eqs. Ijl4|l 
and ljl6JI . Its evaluation yields^ 7 - 

dl 1 f 

^F = G ^- / du{-df/du) [-nvImT ls (oj + eV)} 

(18) 

for the differential conductance. Here /(w) is the Fermi 
function (oj is the energy measured from the Fermi level) , 

and Ti s is the t-matrix for the particles of channel 7 = 1 
[evaluated with the equilibrium Hamiltonians Eq. or 
Eq. JTD}]. The t-matrix is related to the exact retarded 
Green function Gks,k's' = 5 ss >Gks,k's of these particles 
according to 

Gks,k's = G° + G°T ls G°,, G° = (u - & + My 1 . 

Here we took into account the conservation of the total 
spin, which implies that Gk s ,k's' is diagonal in s, s'. In 
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our model with t\ independent of k (and, consequently, 
Ji independent of k and k'), the t-matrix is also inde- 
pendent of fc, k 1 . Note that the linear response (V — > 0) 
counterpart of Eq. (|18fl , the linear conductance 



(20) 



remains valicUi for an arbitrary relation between tz and 
tR, in which case Go = (2e 2 /h) sin 2 (2# )- 



IV. TRANSPORT AT FINITE TEMPERATURE 
AND BIAS 

Equation (|18|) provides a direct link between the mea- 
surable quantity, the differential conductance dl /dV, and 
the properties of the 2CK model, Eq. (|10|l . In the 
channel-symmetric case J\ = J2 = J the NFL behav- 
ior manifests itself in a nonanalytic dependence of the 
t-matrix on energy and temperature^, which leads to a 
rather unusual scaling of the differential conductance at 
low bias and temperature (|eV|,T <C Tk)' 



1 dl _ 1 

C^dV ~ 2 




2CK 



\eV[ 
ttT 



(21) 



The function faCK (%) here is a universal (parameter-free) 
scaling function 9 with the asymptotes 



F 2 CK (x) 



1 + ex 2 

3 r- 



X < 1 

a; > 1 



(22) 



where c is a numerical coefficient of the order of 1. The 
limit eV/T -> of Eq. (J2U> yields 



G = G - 



1 



1 - 



(23) 



for the linear conductance (this result is valid for arbi- 
trary value of thltji). The estimate^ of the Kondo tem- 
perature Tk introduced in Eqs. I|21|) and (|23[) readsiS 

T K ~ Eo(yJ)e^l J \ E = mm{E 1 ,E 2 }. (24) 

The validity of Eqs. (Jjnj and J23Jl is limited by the 
requirements that both the Zeeman energy B and the 
level spacing SE2 are small compared to T, and that the 
exchange constants in Eq. (|10|l are equal to each other: 
Ji = J2. When the system is tuned away from this spe- 
cial point, at a finite 



A = v Ji — VJ2 



(25) 



the conductance changes drastically. In the ideal case 
of T = and 8E2 = 0, the conductance has a step-like 
dependence on A, 



G(A) = G„0(A). 



(26) 



The discontinuity in Eq. (|26|l reflects a quantum phase 
transition between two different Fermi liquid (FL) states, 
in which the spin of the dot 1 forms a singlet with ei- 
ther the collective spin of the electrons in the leads (FL1, 
A > 0) or with that of the dot 2 (FL2, A < 0). At the 
critical point A = 0, the system exhibits NFL behavior 
down to T = 0. In agreement with the general theory 
of quantum phase transitions^, the T — > asymptotics 
at |A| corresponds to the FL, whereas the NFL be- 
havior Eq. \2'M) is preserved at temperatures well above 
certain A-dependent crossover scale Ta, see Fig- EI By 
the same token, the step in the A-dependence of G(A), 
Eq. H26|) . is smeared at finite temperatures. 




FIG. 2: Quantum phase transition between two FL states. 
The NFL behavior is preserved at |A| 7^ 0, provided the tem- 
perature exceeds the crossover scale Ta, see Eq. 13H . The 
width At of step in the conductance G(A) scales with tem- 
perature as vT, see Eq. (1341 . 

In order to estimated the energy scale Ta we consider 
the renormalization group (RG) flow of the effective ex- 
change constants as the high-energy cutoff D is reduced 
from its initial value D ~ Eq. We are interested in the 
case when the bare value of A is small, 



where 



|A| « J, 



J = u{,h + J 2 )/2 ■ 



(27) 



The evolution of the effective coupling constants J* . A* 
with the decrease of D is then described by the Poor 
Man's scaling equations^ 



dj* 



*\2 



dA* 



d( ^ ' ' d( 
with the initial conditions 



(28) 



J*(D )=J, A*(A,)=A. 

Equations l(25)l are valid as long as A* <C J* <C 1 and 
yield the relation A*/A = (J* /J) 2 . By the time J* has 
grown to be of the order of 1 at D ~ Tk, the value of A* 
characterizing the channel asymmetry reaches 



A*(T K ) ~ A/ J 2 



(29) 



This can be viewed as the initial (at D ~ Tk) value of the 
coupling constant of the relevant -> 2 ® channel-symmetry- 
breaking perturbation. The perturbation will eventu- 
ally drive the system away from the 2CK fixed point at 
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D — > 0. However, if A*(Tr-) <C 1, then one expects the 
behavior of the system in a broad range of energies to be 
still governed by the vicinity of the 2CK fixed point. The 
channel anisotropy is a relevant operator with scaling di- 
mension 1/2, see Ref. 12(1 Hence, the dependence of the 
corresponding coupling constant A* on D is described by 



A*{T K ) 



Tk 
D 



1/2 



(30) 



The condition A*(T A ) 
gives the estimate 



1, together with Eq. ffity. then 



T A ~ [A*(T K )] 2 T K ~ (A 2 / J 4 ) T K . 



(31) 



The RG flow stops at D ~ max{T, |eV|}. Conse- 
quently, at max{TA, \eV\] <C T <C Tr-, the channel 
asymmetry yields a small correction to the conductance 
Eq. (|23[1 . The correction is first order in the corre- 
sponding perturbation, hence proportional to A*(T) ~ 
(T A /T) 1 /2 ; an d its sign is determined by the sign of A: 



5G/G <x sgn(A) 



Ta 
T 



1/2 



(32) 



On the other hand, for T, \eV\ <C Ta the system is a 
Fermi liquid, see Fig. Substitution of the t-matrix in 
the form 



-vwlm T ls = 0(A) - sgn(A) 
[cf. Ref. 13 into Eq. O then yields 



3co 2 + tt 2 T 2 



2T 2 



1 dl 

G^dV 



6(A) - sgn(A) 



Ta) [ 2 \nT y 

'(33) 

Again, the linear response (V — > 0) counterpart of 
Eq. (|33|l is valid at any ratio thjtR- The temperature 
dependence of the linear conductance at fixed small val- 
ues of A is sketched in Fig. [3] 

According to Eq. 133JI . corrections to the zero- 
temperature limit of the linear conductance, the step- 
function Eq. I|26l) . are quadratic in temperature - a typ- 
ical Fermi-liquid result. At a finite temperature, the 
step-function is smeared, see Fig. [21 The characteristic 
width At of the smeared step at temperature T is es- 
timated by solving the equation Ta ~ T for A, which 
results in 



A T ~ J 2 VT/Tk. 



(34) 



This "sharpening" of the A-dependence of the linear con- 
ductance with decreasing temperatures (see Fig. [2J can 
be regarded as a "smoking gun" for non-Fermi-liquid be- 
havior. In fact, it might be easiest to first identify un- 
ambiguously the step-like dependence of the conductance 
on A and then use it to tune the device precisely to the 
symmetry point in order to observe the distinctive scaling 



G/G 




FIG. 3: Sketch of the temperature dependence of the linear 
conductance at fixed values of A and Tk. For A < the de- 
pendence is nonmonotonic, with a maximum at T ~ VTWTk. 
At T > T K the conductance scales as G/Go oc [ln(T/TA')]~ 2 , 
see, e.g., Ref.ll7l 



of the differential conductance Eq. I|21|) . Experimentally, 
the value of A is controlled^ by the asymmetry of the 
conductances of the corresponding tunneling junctions, 
which in turn are controlled by the potentials V g on the 
gates forming the junctions. In the vicinity of the sym- 
metry point, the dependence of G on V g should have the 
form of a smeared step-function, whose width SV g should 
scale with temperature as VT, see Fig. [5] 



V. LINEAR CONDUCTANCE AT A FINITE 
MAGNETIC FIELD 

The magnetic field dependence of the linear conduc- 
tance across the device also reveals the critical behavior. 
In this Section we study the dependence G(B) at T = 
in the vicinity of the quantum critical point A = 0. We 
consider only the Zeeman effect of the magnetic field, 
and dispense with its orbital effect (this is an adequate 
approximation for a field applied in the plane of a lateral 
quantum dot device). 

Similar to the effect of a finite temperature, see Fig- El 
the application of a magnetic field at small A results in 
a crossover from the limiting FL behavior at B — * to 
NFL intermediate regime at higher fields B > B&- As 
before, the crossover scale B& can be estimated^ from 
RG arguments. The scaling dimension 2 - of the operator 
S z in Eq. (jTU|l at the 2CK fixed point is 1/2. Accordingly, 
when the high energy cutoff D is lowered, the effective 
splitting of the impurity levels B* evolves according to 



B*{D)/D 
B*(T K )/T K 



1/2 



(35) 



with the initial condition B*(T K ) - B. The RG flow 
Eq. I|35l) terminates once B* has grown to become of the 
order of D, or when D reaches the value Ta, whichever 
occurs at a higher value of D. The first of the two con- 
ditions corresponds to the limitation on the NFL behav- 
ior set by the Zeeman splitting, while the second one is 
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due to the channel anisotropy. Therefore, the crossover 
scale B A can be estimated as that field B ~ B*(Tk) in 
Eq. JSU, at which B*(D) ~ D and D ~ T A simulta- 
neously. Using Eqs. (|35|) and (|31[) . we find the relation 
between the crossover fielo— , the crossover temperature 
Ta, and the channel anisotropy parameter A 



B A ~ VTaTk- (\A\/J 2 )T k . 



(36) 



Note the difference between the A-dependence of the 
crossover temperature Ta [Eq. Q31[)] and the crossover 
field B A . 

Having found the crossover scale B A , next we investi- 
gate the dependence of the conductance G on the field 
B. First of all, we note that at A ^ the low-energy 
properties of the Hamiltonian Eq. (|10() are those of a 
Fermi liquid^ The effect of any local perturbation, such 
as the exchange interaction with the spin of the dot 1 in 
Eq. (|10|l , on the ground state of the Fermi liquid is com- 
pletely characterized by the scattering phase shifts S 1S at 
the Fermi level. (Recall that s = ±1 for spin- up/down 
and 7 = 1,2 labels the two channels.) The t-matrix that 
enters Eq. I|20(l is then given by the standard scattering 
theory expression 



™T 7S (0) = - (e 



2i8„ 



2i 



1) 



(37) 



Obviously, the phase shifts are defined only mod tt (that 
is, 8~f S is equivalent to i5 7S +7r ). The ambiguity is removed 
by setting the values of the phase shifts corresponding to 
J 7 = in Eq. I|10l) to zero. With this convention, the 
invariance of the Hamiltonian (|10fl with respect to the 



particle-hole transformation t/j-yks — > ij) _ h 
into the relation 



5 1S + <5 7i _ s = 



translates 



(38) 



for the phase shifts, which suggests a representation 



(39) 



Substitution of Eqs. and (jSHJ into Eq. J20J) yields 
G/G = lj2 sin2 Sls = sin2 Sl ( 4 °) 



for the linear conductance at T = 0. In the limit 
B/Tk — * +0 and at A ^ 0, the ground state of the 
Hamiltonian i|l(J|) is a singlet. Therefore, the total spin 
in a very large but finite region of space surrounding the 
dot 1 is zero. By the Friedel sum rule, this implies re- 
lation sS ls = 7r. Taking, in addition, Eq. 139|) into 
account, one obtains relation 



Si 



tt/2, 



(41) 



valid at any value of B/B A , as long as B <C Tk- 

Below the crossover, B <C B A , the values of the phase 
shifts are determined by the vicinity of the stable Fermi- 
liquid fixed points^, Si = tt/2, 5 2 — at A > and 



5\ = 0, S2 — tt/2 at A < 0. Substitution of these values 
into Eq. (|4U[) then yields Eq. I)26|) for the conductance. 
The corrections to the fixed point values of the phase 
shifts are linear in B/B A , 

S 1 = tt/2 -S 2 = (tt/2)0(A) - sgn(A)(B/fl A ), (42) 
yielding 

G/G = 9(A)-sgn(A)(B/B A ) 2 , B « B A (43) 
[cf. Eq. Pty]. 

Above the crossover, i.e., for B A <C B <C Tk, the 
departure of the phase shifts from the 2CK fixed point 
values (5i 2 = tt/4 is controlled by the properties of the 
fixed point. To account for a finite value of B/Tk, we 
generalize Eq. (|41|l : 

S 1 +S 2 =tt [1/2 + M(B)] . 

The zero-temperature magnetization AI(B) here is 
known exactly from the Bethe-ansatz solutionSi^i. Us- 
ing the asymptoteSi M(B) cx (B/Tk) \n(T K /B), we find 

*i = ~+asgn(A)^-&Aln^. (44) 

Here a and b are positive numerical coefficients of the 
order of 1. The second term on the right-hand side 
of Eq. (|42(l is the first-order correction in the channel- 
symmetry-breaking perturbation. This correction is sim- 
ilar to Eq. Q32f) with temperature T replaced by the en- 
ergy scale D*(B) - B 2 /T K at which the RG flow defined 
by Eq. terminates. Equations (glj and (gUJ> yield the 
asymptote of the conductance at B A <B< Tk, 

G 1 , . , B A B Tk , . _ n 

— = -+ asgn(A)— -b— In— . (45) 

The shape of G(B) is qualitatively similar to that of 
G(T), see Eqs. I^,©, and 63, although the precise 
functional form is rather different. 

Interestingly, in the case of small channel anisotropy, 
Ta <C Tk, there is an approximate symmetry with re- 
spect to the change of sign of A: 



G(B, A) + G(B, -A) = 2G(B, A — > 0). 



(46) 



Note that this relation is valid at any B/Tk, provided 
that T A /T K < 1. 

Strictly speaking, the consideration of this Section is 
applicable only at zero temperature. However, the results 
Eqs. H43[) and l|45|l remain validi^ as long as 



T < B 2 /T K 



(47) 



At higher temperatures the conductance is described by 
the corresponding expressions of Section Hvl As follows 
from Eqs. (|23|l and l|45|l . the limiting value of the lin- 
ear conductance at the 2CK fixed point, G = Go/2, is 
independent of the order in which the limits B — > 0, 
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T — > are taken2£. Hence, the crossover between the 
field-dominated regime, see Eqs. l|4*5|) and and the 

temperature-dominated one, see Eqs. (J23J, and (fTTT^fl . 
is expected to be smooth and featureless. 

For arbitrary values of T A /Tx, the detailed magnetic 
field dependence of the phase shifts at the Fermi level 
can be studied using the numerical renormalization group 
(NRG) 24 . In this approach one defines a sequence of 
discretized Hamiltonians and diagonalizes them itera- 
tively to obtain the finite-size spectrum of the model. 
In the Fermi liquid case (A ^ 0) knowledge of the finite- 
size spectrum is sufficient to identify unambiguously the 
phase shifts^S. 
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FIG. 5: Dependences of the crossover scales Ba and Ta on 
the asymmetry parameter A = Aj J 2 . 
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FIG. 4: The phase shifts for the 2CK model at different values 
of the channel asymmetry parameter A — A/J 2 . The upper 
(lower) curves represent 81(62). 

In Fig.QJ we plotted the phase shifts #1,2 as a function 
of B for different values of the parameter A = A/ 'J 2 > 
that characterizes the asymmetry between the channels. 
We estimate the crossover scales 1 ^ Tk and B A as the two 
values of B in Fig. 0] at which the phase shift S 2 equals 
tt/8. In order to verify the relation B A /Tx ~ A, see 
Eq. (|36[) . we plotted Ba vs A on the left panel in Fig. [3] 
The NRG data also allow us to estimate the scale Ta, 
see Eq. (ETTJ) . as the energy scale at which the first excited 
state of the NRG spectrum has reached the halfway mark 
of its crossover evolution between the corresponding two 
fixed point values, see Fig. [5] right panel. The NRG data 
are very well described by B A /T K 0.5A, T A /T K w 
4A 2 , in agreement with Eqs. (|3fcil> and above. 

Having extracted the phase shifts, we are able to cal- 
culate the linear conductance from Eqs. (|40J) and (|46ll . 
see Fig. [S] As expected, the conductance develops a sig- 
nature of a plateau at intermediate values of the field 
B A < B < T K . At very high fields, B > T K , the con- 
ductance scales with B as 1/ \t?(B /Tx). 

As usual in NRG calculations, we measured all ener- 
gies in units of the bandwith D. In order to avoid the 
disturbing finite bandwidth effects, we used two different 
coupling constants for the high- and low-field regimes: 
one set of data, that includes the B 3> B A regime, was 
obtained using J — 0.075, while another set of data, 



FIG. 6: Field dependence of the conductance at different val- 
ues of the asymmetry parameter A — A) J 2 . The upper 
(lower) curves correspond to A > (A < 0). 



which includes the B -C Tk regime, was obtained using 
J = 0.15. The two sets were combined by rescaling the 
magnetic field in units of the Kondo temperature, result- 
ing in a set of continuous curves, as shown in the figures. 
The overlap of the two sets of data at intermediate fields 
confirms that in this regime the accuracy of our numer- 
ics is remarkably good. Based on the dependence on the 
finite system size, we estimate the relative error of the 
calculated phase shifts to be of the order of 2%. (The 
worst case is the low field part of the A = curve, be- 
cause of the extremely fragile nature of the intermediate 
NFL fixed point.) 



VI. EFFECT OF POTENTIAL SCATTERING 

So far we concentrated on the particle-hole symmetric 
model. In general, however, such symmetry is absent. It 
is violated by the presence of higher energy levels in dot 
1, and also by deviations of the dimensionless gate volt- 
age N from an integer value. In the absence of particle- 
hole symmetry, the effective Hamiltonian <|10[) acquires 
additional terms leading to potential scattering. Taking 
into account that the interchannel scattering is blocked 
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at energies well below -Ei,2, we can write this additional 
perturbation as 

7=1,2 kk's 

Including H p into our considerations leads to a modifi- 
cation of the limiting values of the conductance in the 
Fermi-liquid and 2CK fixed points. The dependences of 
dl/dV on A, V, T and B, however, remain the same apart 
from acquiring a constant background contribution G e i 
due to elastic cotunneling. Here we illustrate this for 
a specific example of the zero-temperature magnetocon- 
ductance. 

The potential scattering yields finite spin-independent 
phase shifts <5? = — arctan(7rz/V~ 7 ) even if J 7 in Eq. (|1(J|> 
are set to 0. This can be accounted for by a proper 
modification 25 of Eq. lj3"9"jl . 

5 7S = 5° + s<5 7 , (49) 

where the dependence of 5 7 on B and A is described by 
the "particle- hole symmetric" expressions (|42[1 and (14411 . 
Substitution of the phase shifts in the form of Eq. 1491) 
into Eq. ij4"U|> results ini& 

G(B, A) = G el + G Q F[B/B A , B/T K , sgn(A)], (50) 

where G e i — Go sin 2 5®, the function F is a universal 
function with asymptotes given in Eqs. I|43[l and 1)45(1 . 
and Go — Go — 2G e i. Note that the limiting value of 
the conductance at the 2CK fixed point, G e i + Go/2, lies 
precisely half-way between the two Fermi-liquid limits, 
G e i and G e i + Go, and that Eq. I|46|l remains valid even 
in the presence of the potential scattering Eq. 

VII. DISCUSSION 

The low-tempcrature properties of a quantum dot de- 
vice normally are well described by Fermi liquid the- 
ory. The special two-dot structure proposed in Ref. |l| 
allows, however, for NFL behavior at a special point in 
the space of parameters of the device. In the context of 



the physics of quantum phase transitions, this point can 
be viewed as a critical point separating two Fermi liq- 
uid states. In this paper, we developed a detailed theory 
of the transport properties near such a quantum criti- 
cal point. Our theory offers a strategy for tuning the 
device parameters to the critical point characterized by 
the two-channel Kondo effect physics, by monitoring the 
temperature dependence of the linear conductance, see 
Sec. IIVI Further confirmation of the 2CK behavior may 
come from the measurements of the differential conduc- 
tance, which must display universal behavior, see Sec. IIVI 
We also investigated the effect of magnetic field and of 
potential scattering on the conductance in the vicinity 
of the quantum critical point, see Sees. IVland lVll The 
Zeeman splitting allows one to investigate the finite-field 
crossover between the Fermi liquid and NFL behavior of 
the conductance. In the vicinity of the NFL point, the 
linear conductance of the device depends on the magnetic 
field and temperature only via two dimcnsionlcss param- 
eters, T/Ta and B/Ba; the dependence of Ta and B\ on 
the channel asymmetry A is given in Eqs. (|31|) and 13011 . 
Note also that potential scattering does not destroy the 
2CK behavior, but merely renormalizes the magnitude 
of the Kondo contribution to the conductance. A finite 
level spacing in the larger dot 8E2, however, is a hazard. 
At temperatures below 6E2 the two-dot device inevitably 
enters into the conventional Fermi-liquid regime. 
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